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ABSTRACT 

We have undertaken the largest systematic study of the high-mass stellar initial mass function (IME) to date 
using the optical color-magnitude diagrams (CMDs) of 85 resolved, young (4 Myr < t < 25 Myr), intermediate 
mass star clusters (10^-10"^ Mq), observed as part of the Panchromatic Hubble Andromeda Treasury (PHAT) 
program. We fit each cluster’s CMD to measure its mass function (ME) slope for stars > 2 Mq. By modeling 
the ensemble of clusters, we find the distribution of ME slopes is best described by F = +1.45 ;^q with a very 
small intrinsic scatter. This model allows the MF slope to depend on cluster mass, size, and age, but the data 
imply no significant dependencies within this regime of cluster properties. The lack of an age dependence sug¬ 
gests that the MF slope has not significantly evolved over the first ~ 25 Myr, and provides direct observational 
evidence that the measured MF represents the IMF. Taken together, this analysis — based on an unprecedented 
large sample of young clusters, homogeneously constructed CMDs, well-defined selection criteria, and consis¬ 
tent principled modeling — implies that the high-mass IMF slope in M31 clusters is universal. The IMF has 
a slope (r =+1.45!°:“) that is slightly steeper than the canonical Kroupa (+1.30) and Salpeter (+1.35) values, 
with no drastic outliers in this sample of nearly 100 clusters. Using our inference model on select Milky Way 
(MW) and EMC high-mass IMF studies from the literature, we find Fmw ^+1.15±0.1 and Flmc ^ +1.3 ± 0.1, 
both with intrinsic scatter of ^0.3-0.4 dex. Thus, while the high-mass IMF in the Local Group may be univer¬ 
sal, systematics in literature IMF studies preclude any definitive conclusions; homogenous investigations of the 
high-mass IMF in the local universe are needed to overcome this limitation. Consequently, the present study 
represents the most robust measurement of the high-mass IMF slope to date. To facilitate practical use over 
the full stellar mass spectrum, we have grafted the M31 high-mass IMF slope onto widely used sub-solar mass 
Kroupa and Chabrier IMFs. The increased steepness in the M31 high-mass IMF slope implies that commonly 
used UV- and Ha-based star formation rates should be increased by a factor of ^1.3-1.5 and the number of 
stars with masses > 8 Mq are ^25% fewer than expected for a Salpeter/Kroupa IMF. 

Keywords: stars: luminosity function, mass function — Hertzsprung-Russell and C-M diagrams — galaxies: 
star clusters — Local Group — galaxies: star formation 


1. introduction 

The stellar initial mass function (IMF), for M > 1 Mq 
plays a central role in a wide variety of astrophysics. The 

* Based on observations made with the NASA/ESA Hubble Space Tele¬ 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are associated with pro¬ 
gram #12055 

^ Astronomy Department, Box 351580, University of Washington, Seat¬ 
tle, WA, USA; dweisz@uw.edu 

^ Center for Cosmology and Particle Physics, New York, University, 4 
Washington Place, New York, NY 10003, USA 

^ Raytheon Company, 1151 East Hermans Rd, Tucson, AZ 85756, USA 

^ Max Planck Institute for Astronomy, Koenigstuhl 17, 69117 Heidel¬ 
berg, Germany 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, 
Cambridge, MA 02138, USA 

^ Department of Astronomy, University of Michigan, 500 Church Street, 
Ann Arbor, MI 48109, USA 

^ Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, 
MD 21218, USA 

^ Universitat Heidelberg, Zentrum fur Astronomie, Institut fur Theo- 
retische Astrophysik, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany 

^ Department of Astronomy, University of California at Santa Cruz, 
1156 High Street, Santa Cruz, CA, 95064, USA 

Department of Physics and Astronomy, University of Utah,Salt Lake 
City,UT 84112, USA 

Minnesota Institute for Astrophysics, University of Minnesota, 116 
Church Street SE, Minneapolis, MN 55455, USA 
Hubble Fellow 


relative numbers of such stars is essential to interpreting the 
stellar populations of star-forming galaxies across cosmic 
time, testing and validating theories of star formation, con¬ 
straining chemical evolution models, the formation of com¬ 
pact objects and gravitational waves, and understanding the 
interplay between stars and gas (e.g., Schmidt 1959; Tins¬ 
ley 1968; Searle et al. 1973; Talbot & Arnett 1973; Ostriker 
et al. 1974; Audouze & Tinsley 1976; Scalo 1986; Kenni- 
cutt 1998a; Elmegreen 1999; Kroupa 2001; Massey 2003; 
Alves et al. 2007; Eardal et al. 2007; McKee & Ostriker 2007; 
Wilkins et al. 2008; Pfiamm-Altenburg & Kroupa 2009; Ken- 
nicutt & Evans 2012; Narayanan & Dave 2012; Conroy 2013; 
Belczynski et al. 2014; Krumholz 2014; Madau & Dickinson 
2014). 

1.1. Current State of the IMF Above ^ I Mq 

Our best constraints on the high-mass IME come from stud¬ 
ies in the Milky Way (MW) and nearby galaxies. Eor the clos¬ 
est galaxies, it is possible to count individual stars and make a 
direct measurement of the high-mass IME slope. While mea¬ 
suring the IME from direct star counts is straightforward in 
principle, it is far more complicated in practice. Accurate 
measurements are challenging due to a variety of observa¬ 
tional and physical effects including completeness, dynami¬ 
cal evolution, distance/extinction ambiguity, stellar multiplic¬ 
ity, the accuracy of stellar evolution models, degeneracies be- 
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tween star formation history (SFH) and the IMF slope, and a 
paucity of massive stars in the nearby universe, all of which 
complicate the translation of the present day mass function 
(MF) into the IMF (e.g., Lequeux 1979; Miller & Scalo 1979; 
Scalo 1986; Mateo 1993; Massey et al. 1995a; de La Fuente 
Marcos 1997; Maiz Apellaniz & Ubeda 2005; Elmegreen & 
Scalo 2006; Maiz Apellaniz 2008; Ascenso et al. 2009; Krui- 
jssen 2009; De March! et al. 2010; Portegies Zwart et al. 2010; 
Bastian et al. 2010). 

Observations of more distant, unresolved galaxies alleviate 
some of these challenges. For example, by studying an entire 
galaxy, effects of stellar dynamics (e.g., migration, mass seg¬ 
regation) and distance/extinction confusion are less important. 
Further, distant galaxies span a more diverse range of environ¬ 
ments than what is available in the very local universe, offer¬ 
ing better opportunities to measure the high-mass IMF as a 
function of environment (metallicity, mass, redshift, etc.; e.g., 
Baldry & Glazebrook 2003; Fardal et al. 2007; Dave 2011; 
Narayanan & Dave 2013). However, the unresolved nature 
of these observations severely limits the precision and accu¬ 
racy of a high-mass IMF measurement. Most integrated ob¬ 
servations have insufficient leverage to mitigate the SFH-IMF 
degeneracy, and are forced to make simplified assumptions 
about a galaxy’s recent SFH (e.g., constant, single short burst) 
and dust (e.g., a single value for an entire galaxy) to provide 
any constraint on the high-mass IMF (e.g.. Miller & Scalo 
1979; Elmegreen & Scalo 2006). Although studies of unre¬ 
solved clusters may provide a new avenue for IMF studies 
(e.g., Calzetti et al. 2010; Andrews et al. 2013), high-mass 
IMF studies using integrated light typically provide coarse 
characterizations rather than precise measurements. 

As a result of these challenges, our current knowledge of 
the high-mass IMF remains remarkably poor. The widespread 
adoption of a ‘universal’ Salpeter (F = -hi.35; Salpeter 1955) 
or Kroupa (F = -hl.3; Kroupa 2001) IMF above 1 Mq, pre¬ 
sumed not to vary with environment has given the impression 
that the high-mass IMF is a solved problem (e.g., Kennicutt 
1998a). However, compilations of literature IMF studies from 
resolved star counts indicate an “average value” of F = 1.7 
with a scatter >0.5 dex in IMF slope measurements that do 
not clearly support a ‘universal’ IMF (e.g., Scalo 1986, 1998; 
Kroupa 2002; Scalo 2005)^"^. There also does not appear to be 
any unambiguous systematic dependence of IMF slope on en¬ 
vironment and/or initial star-forming conditions (e.g., Bastian 
et al. 2010), which has significant implications for how stars 
form in a wide variety of physical conditions (e.g., Elmegreen 
1999; McKee & Ostriker 2007; Hopkins 2013; Krumholz 
2014). 

Controversy over the form and universality of the IME have 
persisted since the first studies in the 1950s. Eollowing the 
seminal work of Salpeter (1955), a handful of studies re¬ 
ported similar IMEs in a number of MW clusters and field 
environments (Jaschek & Jaschek 1957; Sandage 1957; van 
den Bergh 1957), prompting early speculation about the uni¬ 
versality of the “original mass function” (e.g.. Limber 1960; 
van den Bergh & Sher 1960). However, this early agreement 
with Salpeter gave way to large scatter in IME measurements 
beginning in the 1970s, as several studies with more sophis¬ 
ticated completeness corrections and updated stellar physics 

Although a ‘Kroupa IMF’ has a high-mass slope of F = 1.3, Kroupa 
(2002) argue that the slope should probably be close F ~ 1.7, as most histor¬ 
ical IMF studies do not account for unresolved binary stars, which, in princi¬ 
ple, cause measured MF slopes to be shallower than true IMF slopes. 


highlighted the shortcomings in the first generation of IME 
studies (e.g., Taff 1974; Lequeux 1979; Tarrab 1982). Eur- 
thermore. Miller & Scalo (1979) provided a comprehensive 
reassessment of the MW field star IME, in light of the previ¬ 
ously unaccounted for SEH-IME degeneracy. They suggested 
a high-mass IME of the with F = 1.5 for stars 1 <M/M0<1O 
and F = 2.3 for stars M/Mq>10. However, Kennicutt (1983) 
showed that such a steep IMF slope above ^10 Mq was in¬ 
compatible with observations of B-V colors and Ha equiv¬ 
alent widths in nearby disk galaxies, suggesting that extrapo¬ 
lating the value of F = 1.5 from 1-100 Mq provided a better 
match to the data. This conclusion was further solidified in 
Kennicutt et al. (1994). 

The number of observational high-mass IMF studies 
peaked between the late 1980s and early 2000s. This golden 
age followed the comprehensive IMF review by Scalo (1986), 
which concluded that systematic uncertainties were the dom¬ 
inant source of confusion in our poor knowledge of the IME 
Shortly after, the proliferation of CCDs enabled significant 
progress in reducing these systematics, primarily by allow¬ 
ing for quantitative completeness corrections via artificial star 
tests (e.g.. Stetson 1987; Mateo 1993). The 1990s saw a fiurry 
of qualitatively new IME studies including star counts in the 
LMC, spectroscopy of starburst galaxies, and IME constraints 
from stellar abundance patterns (see Leitherer 1998; Massey 
1998; Gilmore 2001 for extensive discussion of the literature 
in these areas). The sheer number and richness of IME stud¬ 
ies from this era precludes a detailed history in this paper, and 
we instead refer readers to a number of excellent papers that 
chronicle the history of the IME (e.g., Scalo 1986; Massey 
1998; Kennicutt 1998a; Elmegreen 1999; Eisenhauer 2001; 
Gilmore 2001; Massey 2003; Scalo 2005; Elmegreen 2006; 
Zinnecker 2005; Bastian et al. 2010). 

To briefly summarize, many studies from this era concluded 
that the high-mass IME in the local universe was not drasti¬ 
cally different from Salpeter over a wide range of environ¬ 
ments (e.g., Massey 1998; Leitherer 1998; Gilmore 2001; 
Kroupa 2001; Massey 2003). However, there was no clear ex¬ 
planation for the observed scatter in the measured IME slopes, 
which contradicts a truly universal high-mass IME. Some have 
argued that measurement uncertainty and systematics are the 
dominant source of scatter (e.g., Elmegreen 1999; Kroupa 
2001, 2002; Massey 2003). However, without systematically 
re-visiting each study, it is challenging to dismiss the scat¬ 
ter as solely a measurement uncertainty. As a clear counter 
example, there are well-known studies of star clusters with 
some of the same physical properties (e.g., age, mass) that 
have been analyzed identically, but have very different high- 
mass IME slopes (e.g., Phelps & Janes 1993), reinforcing that 
some of the scatter may be physical in nature. At the same 
time, on the whole, >75% of resolved star IME studies have 
significantly underreported uncertainties relative to the theo¬ 
retically possible precision (Weisz et al. 2013), which under¬ 
lines the difficulty in disentangling systematics from intrinsic 
IME variation using literature IME studies. 

Since the early 2000s, studies of the high-mass IME have 
diminished in number, as attention turned to the sub-solar 
regime (e.g., Chabrier 2003; Scalo 2005; Zinnecker 2005). 
However, focus on the high-mass IME was recently rekindled 
following several reports of a systematic deficit of massive 
stars in lower-luminosity galaxies, based on observations of 
integrated broadband and/or Ha-to-UV fiux ratios (e.g., Hov- 
ersten & Glazebrook 2008; Boselli et al. 2009; Lee et al. 2009; 
Meurer et al. 2009). Yet, whether these observations are in- 
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dicative of a systematically varying IMF (e.g., Weidner & 
Kroupa 2005 ; Pflamm-Altenburg & Kroupa 2009) remains an 
open question, as other mechanisms such as stochastic sam¬ 
pling of the cluster mass function and/or bursty SFHs can also 
explain these observed properties (e.g., Fumagalli et al. 2011; 
da Silva et al. 2012; Weisz et al. 2012). 

Thus, at present, the high-mass IMF is characterized by 
common usage of a universal Salpeter/Kroupa IMF, but also 
boasts a vast body of literature studies which does not unam¬ 
biguously support this as the de facto standard. After decades 
of study, the main limitation in our knowledge of the high- 
mass IMF remains our reliance on a rich, but heterogenous 
set of literature IMF studies, instead of a large, homogenous 
survey of individually resolved young, massive stars. 

1.2. The IMF > 2-3 M© in M31 

A primary goal of the Panchromatic Hubble Andromeda 
Treasury program (PHAT; Dalcanton et al. 2012) is a com¬ 
prehensive exploration of the stellar IMF above ^ 2-3 M©. 
This 828-orbit HST multi-cycle program acquired near-UV 
through the near-IR imaging of ^120 million resolved stars 
in ^ 1/3 of M31’s star-forming disk (Williams et al. 2014), 
providing a dataset of unprecedented size and homogeneity 
for measuring the high-mass IMF. 

The focus of the present IMF study is on M3rs young star 
cluster population. Star clusters provide a more direct path to 
measuring the IMF compared to field populations, which suf¬ 
fer from IMF-SFH degeneracies for stars more luminous than 
the oldest main sequence turnoff (e.g.. Miller & Scalo 1979; 
Elmegreen & Scalo 2006), Further, the differential extinction 
found in the disk of M31 (e.g., Draine et al. 2013; Dalcan¬ 
ton in prep.), and similarly sized star-forming galaxies, adds 
a significant layer of complexity to modeling the field stel¬ 
lar populations. In contrast, star clusters provide a relatively 
simple way to probe the IMF. Their co-eval nature genuinely 
mitigates the effects of SFH, while their vulnerability to fore¬ 
ground differential extinction is minimal due to their small 
sizes. 

In this paper, we focus on determining the MF slopes for 
PHAT clusters that are < 25 Myr old and ^10^ M© in mass 
using resolved star color-magnitude diagrams (CMDs). For 
these clusters, the detectable portion of the main sequence is 
well-populated from ^ 2M© to > 25M© and they are of or¬ 
der ^0.1 relaxation times old, which mitigates the strongest 
effects of mass segregation, making them an ideal sample 
with which to study the high-mass IMF. We model the optical 
CMD of each cluster k to constrain its power-law slope Tk of 
the present day high-mass MF, p(r^), marginalized over other 
parameters of that fit (e.g., cluster age, extinction, binarity). 

Subsequently, we propose a simple model for the ensemble 
properties of the cluster MF, with a mean, F, and a scatter, crp, 
and fit it to the MF probability distribution functions (PDFs) 
of the individual clusters. We also use this model to search for 
possible dependencies of F on cluster mass, age, and size. We 
verify this approach with an extensive set of artificial clusters 
that were inserted into PHAT images, processed, and analyzed 
identically to the real clusters. 

This paper is organized as follows. In §2, we define our 
cluster sample and briefly describe the PHAT photometry, ar¬ 
tificial star tests, and artificial cluster sample. We outline our 
method of modeling the cluster CMDs in §3 and present the 
results for the ensembles of real and artificial clusters in §4. 
Finally, in §5, we discuss our findings in the context of current 
and historical views on the high-mass IMF and illustrate the 


implications of our results 

2. DATA 

2.1. Cluster Sample 

Clusters in PHAT were identified as part of the Andromeda 
Project^^, a collaborative effort between PHAT and the Zooni- 
verse^^ citizen science platform. As described in Johnson 
et al. (2015), cluster candidates were assigned probabilities 
of actually being a cluster based on citizen scientist ranking. 
These probabilities were verified using synthetic cluster re¬ 
sults and comparison with the preliminary catalog, which was 
assembled by professional astronomers (Johnson et al. 2012). 
Extensive details of the cluster identification methodology can 
be found in Johnson et al. (2015). 

For the present study, we have limited our analysis to the 85 
clusters with CMD-based ages < 25 Myr and masses ^10^ 
M© (Beerman in prep.), and half-light radii ^ 2 pc (Johnson 
et al. 2015). We focus on these young, more massive clusters 
in PHAT to avoid the very sparse sampling in the low-mass 
clusters; to limit the impact of stellar crowding in the densest 
clusters; to facilitate the extrapolation of the observable MF 
to an IMF; and to remain in the regime of high completeness 
(>90% for our sample; see Figure 10 in Johnson et al. 2015). 
We will explore each of these issues in detail in our compre¬ 
hensive study of PHAT cluster MFs in a future paper. 

2.2. Cluster Photometry 

Photometry of PHAT clusters was performed using 
DOLPHOT^^, a version of HSTPHOT (Dolphin 2000) with 
HST specific modules, following the description in Williams 
et al. (2014)^^. Slight changes to the survey-wide photometry 
include reducing only the optical photometry (F475W, simi¬ 
lar to SDSS g; and F814W, similar to I-band), which, for the 
purposes of this study provide the most information. The UV 
imaging is significantly shallower (^ 2 mag) than the opti¬ 
cal, even along the main sequence, resulting in many fewer 
sources detected. The IR also yields far fewer main sequence 
stars as a result of significant crowding due to the lower angu¬ 
lar resolution of the WFC3/IR camera. An image and optical 
CMD of a typical cluster used in this study are shown in Fig¬ 
ures 1 and 2. 

In principle, the full spectral energy distribution of each star 
can increase the degree of precision to which we know its 
mass, and, in turn provide improved constraints on the IMF. 
However, the UV and IR photometry in PHAT clusters typi¬ 
cally contains fewer than ~ 50% the number of stars detected 
in the optical, reducing their statistical utility for constraining 
the MF slope (cf. Weisz et al. 2013). Further, the UV and 
IR bolometric corrections are not as certain as those in the 
optical, particularly for massive, metal rich stars, which can 
introduce systematics into our analysis. Finally, the modest 
gain (at best) in precision on the IMF slope of each cluster 
must be balanced against the computational cost of process¬ 
ing >10 million additional artificial star tests (ASTs). Thus, 
by only reducing and analyzing the optical data, we retain vir¬ 
tually all statistical leverage on measuring the MF slope, but 
save considerably on computational costs. 

However, analysis of the full stellar SEDs are valuable for 
searching for the most massive stars in clusters. In addition to 

http://www.andromedaproject.org 

https : / /www. zooniverse . org 

http : //americano . dolphinsim. com 

http://archive.stsci.edu/missions/hlsp/phat/ 
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Figure 1. The HST/ACS optical color image of AP 94, a typical cluster ana¬ 
lyzed in this paper. The angular high-resolution of HST allows us to resolve 
the individual stars clusters at the distance of M31. For reference, the image 
is ~50pc on each side. 

the slope, the maximum stellar mass plays a fundamental role 
in defining the high-mass IMF, but is unconstrained by the 
approach we use to measure the MF slope (e.g., Weisz et al. 
2013). An investigation of the most massive stars in the PHAT 
area is the subject on an ongoing study and will be presented 
in a future paper. 

For each cluster, we ran ^50,000 ASTs by adding indi¬ 
vidual stars of known luminosities to PHAT images of each 
cluster. The ASTs were uniformly distributed over each clus¬ 
ter’s CMD and spatially distributed approximately according 
to each cluster’s light profile. Extensive testing showed that 
such a spatial scheme was necessary to accurately character¬ 
ize the completeness and photometric uncertainty profiles of 
the clusters, as a simple uniform spatial scheme results in 
overly optimistic completeness and photometric uncertainty 
profiles. 


2.3. Artificial Clusters 

We use extensive sets of artificial clusters to verify the accu¬ 
racy of each component of our MF determinations, including 
photometry, ASTs, and CMD modeling. The artificial clus¬ 
ters were designed to match the physical properties of our 
real cluster sample and span a range of physical parameters 
(age, mass, size, extinction) that encompass the real clusters. 
The artificial clusters were inserted into images at a range of 
galactocentric radii and environments to capture the full range 
of background densities in PHAT. Overall, the artificial clus¬ 
ter tests confirmed that we can accurately recover the true MF 
from clusters in PHAT imaging, along with all other cluster 
physical properties. 


3. METHODOLOGY 

3.1. Measuring the Stellar Mass Function of a Cluster 

To measure the MF of a cluster we construct synthetic mod¬ 
els of its optical CMD using MATCH, a CMD fitting program 
described in Dolphin (2002). Briefiy, for a given stellar evo¬ 
lution library, binary fraction, IMF, age, metallicity, extinc¬ 
tion, and distance, MATCH constructs a simple stellar popula¬ 
tion (SSP) in the desired filter combinations. This synthetic 
SSP is then convolved with completeness and photometric bi¬ 
ases and uncertainties as determined by the ASTs. The result¬ 
ing synthetic CMD is then linearly combined with a suitable 
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Figure 2. Panel (a) - The optical CMD of AP 94. The grey points are below 
the nominal completeness limit, and were not included in the CMD analysis. 
For reference, we indicate the main sequence turn off masses on the left hand 
side and over-plot the best fitting isochrone for this cluster (see §3.1 and 
Figure 3). Panel (b) - the background CMD constructed by selecting all 
stars in an annulus surrounding the cluster that has lOx times the cluster area. 
The black points in this CMD are used to statistically model background 
contribution to the cluster CMD. 

CMD of foreground/background populations, with the rela¬ 
tive contribution of each controlled by a scale parameter. For 
a given set of these parameters, the synthetic Hess diagram, 
the binned number density of stars in the CMD, is compared 
to the observed CMD using a Poisson likelihood function, 
which is necessary to account for the effects of sparse sam¬ 
pling. The procedure is repeated over a large grid of parame¬ 
ters in order to map all of likelihood space. For this analysis, 
we adopted the Padova stellar evolution models (Girardi et al. 
2002; Marigo et al. 2008; Girardi et al. 2010) and we have 
listed the other model parameters and their ranges in Table 1 . 

Of these parameters, we limited the metallicity to a fairly 
small range near Zq . This was done because (a) the cluster 
CMDs are too sparsely populated to provide robust metallicity 
constraints and (b) the current metallicity of M31’s disk is 
known to be approximately solar, with a very weak gradient 
based on HII region and supergiant abundances (e.g., Venn 
et al. 2000; Trundle et al. 2002; Sanders et al. 2012; Zurita & 
Bresolin 2012). 

Similarly, we adopted a fixed binary fraction of 0.35, where 
the mass ratio of the stars is drawn from a uniform distribu¬ 
tion. Historically, stellar multiplicity has been problematic 
for high-mass IMF determinations from luminosity functions 
(e.g., Maiz Apellaniz & Ubeda 2005; Maiz Apellaniz 2008). 
However, for CMD analysis, it is much less of an issue for 
two reasons. First, the addition of color information provides 
clear separation between the single and binary star sequences, 
minimizing confusion between the two. In contrast, when 
only luminosity information is available, the binary fraction 
and IMF slope are essentially degenerate parameters. Sec¬ 
ond, stellar multiplicity is observationally less important for 
high-mass stars compared to sub-solar mass stars. For stars 
above 1 Mq, the majority of companion stars are likely to be 
significantly less massive and less luminous than the primary 
(cf. Kobulnicky et al. 2014), resulting in minimal change to 
the CMDs. 

Within PHAT, we verified that our choice in binary fraction 
does not change our results. For example, in the case of AP 94 
(Figures 2 and 3), we modeled the CMD with binary fractions 
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Table 1 

Parameters for Measuring the MF of an M31 Star Cluster 


Parameter 

(1) 

Range 

(2) 

Grid Resolution 
(3) 

Notes 

(4) 

r 

(-1.2,5) 

0.2 dex 

MF slope for stars > 2 Mq 

Binary Fraction 

0.35 


Fraction of Stars with a Companion 

Log(Age) 

(6.60, 9.0) 

0.1 dex 

Cluster Age 

Ay 

(0.0, 3.0) 

0.05 dex 

Line of Sight Extinction 

(m-M)o 

24.47 


Adopted Distance Modulus to M31 

log(Z) 

(-0.2, 0.1) 

0.1 dex 

Metallicity 


Note. — Parameters and their ranges used for modeling M31 cluster CMDs. The 
adopted distance modulus is from the tip of the red giant branch measurement made by 
McConnachie et al. (2005). 


of 0.35 and 0.70 and found differences in the resulting IMF 
slope constraints to be <0.02 dex. Tests including different 
binary fractions on artificial clusters yielded similarly small 
changes. 

However, we did find that extreme binary fractions of 0 or 1 
resulted in significantly worse CMD fits. In these two cases, 
the model either entirely populated the single star sequence 
(binary fraction = 0) or only sparsely populated it (binary 
fraction =1), which substantially reduced the resemblance of 
the model CMD to the observed CMD. Thus, for the present 
work, we adopted a fixed binary fraction of 0.35, which ap¬ 
pears to be a reasonable value from LMC cluster studies (e.g., 
Elson et al. 1998). 

To model the background, we construct an empirical CMD 
from an annulus around the cluster that is ten times larger than 
the cluster. The large area is necessary to sample all parts 
of the CMD to accurately constrain the background scaling 
parameter. 

Given a set of cluster model parameters (age, extinction, 
metallicity, and MF slope) the resulting model CMD, and the 
empirically derived CMD for fore-/background sources are 
compared to the observed CMD in a likelihood sense (Dol¬ 
phin 2002). The resulting posterior probabilities, p(cAv,r) 
are illustrated in Figure 3 for a typical cluster in the sample. 
The prior on metallicity tightly brackets solar metallicity and 
the cluster mass emerges as a normalization factor, rather than 
a fit parameter. Finally, the smoothness of the p(t,Av,T) of 
each cluster allows us to interpolate in order to increase the 
resolution of our PDFs. 

3.2. Inferring the Intrinsic Distribution of Mass Function 
Slopes From a Set of Clusters 

One challenge for past MF studies was how to combine 
measurements from individual clusters into a broader state¬ 
ment about the MF slope of ‘clusters in general’. Typically, 
the variance weighted mean is used to compute the “average” 
MF from an ensemble (e.g., Scalo 1986). However, this can 
potentially be a biased estimator (e.g., assumptions of nor¬ 
mally distributed uncertainties), and does not necessarily use 
all available information (e.g., the probability that an object is 
actually a cluster). 

Here, we adopt a framework that allows us to infer popula¬ 
tion wide characteristics about the distribution of MF slopes, 
given a set of noisy measurements. This model follows the 
approach and implementation laid out in Hogg et al. (2010) 
and Foreman-Mackey et al. (2014). We briefly describe the 
model below and we refer the reader to those papers for more 
detail. 

As a simple model for the distribution of cluster MFs, we 
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Figure 3. The joint and marginalized distributions of age, MF slope, and 
extinction for a single young M31 cluster AP 94 that resulted from modeling 
the cluster’s optical CMD, which is shown in Figure 2. The dark points in 
each panel reflect the regions of high probability, and the pixel sizes reflect the 
grid resolution as indicated in Table 1 . We use the marginalized distribution 
for r (shaded histogram) in the subsequent modeling of the distribution of 
MF slopes (§3.2). 

assume that it can be described by a normal distribution, 
A/’(r, cTp), with a mean of F and an intrinsic dispersion of ar- 
We initially assume that their parameters are independent of 
age, mass, and size of the cluster. If we assume the probabil¬ 
ity that each object is a cluster, we can then write down 
a Gaussian mixture model likelihood function for the MF of 
clusters, where we simply assume that the MFs of falsely pre¬ 
sumed clusters have a different J\f{Tfalse ^(^rjaise)- This mix¬ 
ture model mitigates the impact of erroneously identified clus¬ 
ters, without forcing a binary decision on which clusters have 
been identified with sufficient confidence. 

More explicitly, the mixture model can be written as 

p(T, I d, Q,) = Q, exp (-(r-r,)/(24)) 

+ il-Qk) exp {-(rfalse-Tk)/(2(7f,^^)) 

( 1 ) 

where T^ is the marginalized PDF for the MF of the k^^ object 
and is the probability that the k^^ object is a cluster. F and 
cTp are the mean and variance of the Gaussian of interest for 
the clusters and F false and ^false nuisance parameters for 

possible contaminating sources. Simply put, with this model, 
the more probable an object is a cluster the more it contributes 
to the MF parameters of interest. In this model context, a 
universal Kroupa MF would be represented by F = 1.30 and 
CTp 0, and universal Salpeter MF by F = 1.35 and crp -A 0. 

This same framework also readily allows us to explore 
whether the mean MF slope depends on any of the other clus¬ 
ter properties such as age, mass, and size. Specifically, we can 
generalize the model to 

— Mk 4 

F (4 ,Mk,rk) = T + am log --+atlog--\- Or log — (2) 

Me tc rc 

where {4,My^,ry^} are the most likely age, mass, and size 
for the object and {4,Mc,rc} are the mean of our clus¬ 
ter sample in each case. We assume priors that are fiat in 
F, log (Jr, Om^ Ot, Or, and sample these parameter’s PDFs 
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Figure 4. 85 young, intermediate mass clusters from PHAT used to measure 
the high-mass IMF in M31. Both panels show the median of the probability 
distribution for for each cluster’s MF slope vs. cluster mass. Points in the 
top panel are color-coded by their age, while those in the bottom panel by 
the probability of being a cluster. In all cases, point sizes are proportional to 
half-light radius and the uncertainties in F reflect the 68% confidence interval 
on the MF slope. 

with the affine invariant ensemble Markov chain Monte Carlo 
(MCMC) program emcee (Foreman-Mackey et al. 2013)^^. 

4. RESULTS 

The analysis laid out in Section 3 proceeds in two steps, 
the first leading to the MF slope measurements for individual 
clusters, p(rk), the second to the characterization of the intrin¬ 
sic, i.e., error corrected, distribution of mass function slopes 
in M31 clusters. 

The result of the first step, the marginalized PDF, p(Tk), for 
all 85 sample clusters is illustrated in Figure 4, showing the 
MF slope and its uncertainty for each cluster versus its mass. 
Points in the top panel are color-coded by age and those in 
the bottom panels by their probability of being a cluster. In 
both cases, the point size is proportional to the cluster’s half- 
light radius. Clusters with larger half-light radii have lower 
crowding and thus more accurate mass function determina¬ 
tion. From a physical standpoint, the size of a cluster may 
also refiect the environmental conditions in which it was bom. 

While we have used a different methodology here than 
Weisz et al. (2013), the error bars on the MF are similar in the 
two approaches, limited in both cases by the number of well- 
detected cluster members and by the range of stellar masses 
they span. For example, consider AP 94, whose CMD and 
derived physical property PDFs are shown in Figures 2 and 

https://github.com/dfm/emcee 
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Figure 5. The joint and marginalized distributions for the distribution of 
cluster MFs described by Equations 1 and 3. In all panels, except for ar, 
the dashed blue lines represent the 16th, 50th, and 84th percentiles of each 
marginalized distribution. The distribution for crp is highly non-Gaussian, 
making the median a poor point estimate. Instead, we use the mode as a 
point estimate, as it reflects the most probable portion of the PDF, and use 
the 68th percentile to represent an upper limit. For referenc e, the solid red 
lines indicate the value of a Kroupa IMF slope in the Gamma panels, the up¬ 
per uncertainty on cfy for a universal IMF as determined by similar analysis 
of artificial clusters, and a value of zero in the remaining panels. This plots 
shows that the mean IMF slope in M31 clusters is steeper than Kroupa, that 
the recovered scatter is consistent with expectations of a universal IMF, with 
a small tail to larger values, and that there are no significant trends between 
the IMF slope and age, mass, or size of the clusters. 

3. Effectively, there are ^120 stars (after background correc¬ 
tion) between ^ 3 and 11 Mq on the upper MS that are used 
for MF determination. Given these numbers, from Figure 7 
in Weisz et al. (2013), we expect a theoretical uncertainty of 
^±0.4 on its MF slope. In comparison, the l-a uncertainty 
from the marginalized MF PDF in Figure 3 is ±0.45, which is 
very close to theoretical expectation. Uncertainties in cluster 
membership, the exact mass range on the main sequence, etc., 
can explain the small difference in these two values. In gen¬ 
eral, the cluster MF slopes are generally less precise for older 
and lower mass clusters, and follow the expected precision 
relationship derived in Weisz et al. (2013). 

Figure 4 also gives a qualitative impression of the ensem¬ 
ble properties of MF measurements. The majority of clusters 
have MF slopes that are near the Kroupa (1.30) or Salpeter 
(1.35) slopes; ^ 70% of the clusters MF slopes are within ^ 
1-cr of Kroupa and ^92% are within 2-a. Remarkably, in 
this large sample of young clusters we find that none have ex¬ 
tremely steep or fiat MF slopes. At least for the parameter 
space covered by our cluster sample, extreme MF slope out¬ 
liers (> 3 - cr) are quite rare. 

The plots also show the lack of strong correlations between 
cluster MF slope and their physical properties. From visual 
inspection, it is clear that the MF slopes do not show signifi¬ 
cant trends as a function of age, mass, size, or probability of 
being a cluster. 

These impressions from Figure 4 are quantified in the sec¬ 
ond step of our analysis, where we propose and constrain a 
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Table 2 

IMF Parameters for Ensemble of Clusters 


Parameter 

Artificial 

Real 

Real 


Clusters 

Clusters 

Clusters 

(1) 

(2) 

(3) 

(4) 

T7 

1 29+®-®^ 

^•^^-0.04 

1 46+®-®'^ 

^•^”-0.07 

1 45+003 
^•^'’-0.06 

crrj 

'-'•'-'■^-0.00 

'-'•'-''’-0.00 


Qm 

0 01 

'-'•'-'^-0.10 


at 

0.05!S;S 

^•^”-0.09 


Gr 

-0.02!S;|o 

'-'•'-'^-0.11 



Note. — Column (2) - Values of the five parameter model for the ensemble IMF 
from Equations 1 and 3. In the case of or, we quote the mode of its marginalized 
PDF, along with the upper 68% confidence interval, as the distribution is highly non- 
Gaussian. Column (3) - Value of the IMF slope by fixing all other parameters to zero, 
i.e., the assumption that the IMF is ‘universal’. The listed values reflect the median and 
surrounding 68% confidence interval. 


5. DISCUSSION 


Figure 6. A comparison between the expected IMF slope, F, from artificial 
clusters inserted into PH AT images (blue) and the real clusters (red), under 
the assumption that the underlying IMF slope is universal, i.e., (jp —0. The 
solid lines indicate the most probable values, while the shaded regions are the 
68% confidence intervals. The uncertainties on the two measurements begin 
to overlap at a 2-cr level. 

simple model for the cluster’s distribution of MFs (see §3.2), 
and constrain its parameters through the comparison with all 
85 p(r,). 

Figure 5 illustrates the result for the five parameters in our 
model: the mean MF slope, F, and its intrinsic dispersion 
(Jr, and the three coefficients that represent linear trends be¬ 
tween the MF slope and cluster_age, mass, and size. This plot 
shows that the mean MF slope (F = 1 .46 ;^qq 7) is steeper than a 
Kroupa IMF, that the scatter (ar = 0.03 ;^qqq) is consistent with 
expectations from a universal IMF as determined by similar 
analysis of the artificial clusters, and that there are no signifi¬ 
cant trends between MF slope and cluster physical properties. 
The small degree of scatter is particularly interesting, giving 
the visual impression of large variation between single clus¬ 
ters Figure 4. The typical cluster has an uncertainty on its MF 
slope of ^ 0.5 dex, but that ensemble scatter scales roughly as 
1 / v^A/ciusters), indicating that a large and homogenous sample 
of clusters are needed to statistically identify a universal IMF. 
Finally, the limit on at implies that there is no evidence that 
the MF slope changes between 4 and 25 Myr. We list sum¬ 
mary statistics for each distribution in Table 2 along with the 
same statistics for our analysis of the population of artificial 
clusters. 

There are two reasons why the artificial clusters results are 
not an exact delta function. First, the finite grid resolution in 
p(Tk) sets a floor of a ^ 0.03; second, as or is a positive defi¬ 
nite quantity, its most likely estimate from noisy data must be 
significantly positive. However, this effect is quite small, as 
the lower limit on the scatter for the artificial clusters is con¬ 
sistent with the minimum possible value, given the resolution 
of our individual cluster PDFs. 

Given that the scatter is consistent with expectations for 
a universal IMF, we can consider a model in which all the 
clusters have the same underlying MF i.e., a -A 0. We find 
an ensemble-wide value of F = 1.45 ;^qq 5, which is consistent 
with results from the hierarchical model. For comparison, the 
same analysis applied to the artificial clusters yields a value 
of F = 1.30 ;^qq 4. Uncertainties on the two mean values be¬ 
gin to overlap at the ^ 2-cr level. We plot the results of the 
‘universal’ MF model in Figure 6. 


The quality of the PHAT data has enabled important 
progress on three issues: (1) quantifying the extent to which 
an MF measurement reflects an initial MF (IMF) estimate; (2) 
determining accurately what the mean IMF slope is in typical 
solar-metallicity (cluster) populations of a large disk galaxy; 
and (3) assessing the extent to which the IMF is ‘universal’, 
in the sense that it shows very little scatter in the slope, and no 
trends of the slope with cluster properties. The fundamental 
conceptual limitation of the present approach is that it only 
applies to clusters that remain grouped for longer than 4 Myr, 
i.e., they have not (yet) disrupted. Our conclusions are based 
on measurements of the high-mass IMF that combine a prin¬ 
cipled and systematic constraint of each cluster’s MF slope, 
from CMD analysis, with a probabilistic modeling framework 
to analyze the distribution of MF slopes. 

Overall, we find that the high-mass IMF for clusters in M31 
is consistent with a universal IMF. That is, as discussed in §4 
and illustrated in Figures 5 and 6, the favored high-mass IMF 
model has a steeper mean slope than Salpeter/Kroupa with a 
small scatter that is inline with expectations for a universal 
IMF. Importantly, we also show that the mean MF slope in 
clusters of different ages changes by AF <0.1 per decade of 
age.This suggests that the MF slope is unchanged over the first 
^ 25 Myrs, and provides direct observational evidence that 
the measured MF represents the IMF. In this sense, our result 
represents an important reference point for understanding the 
potential evolution of the cluster MF with cluster physical pa¬ 
rameters (age, mass, size) and local star-forming conditions, 
both of which will be explored in Weisz (in prep.). 

Of these results, the most robust determination is the mean 
IMF slope. Even if we assume that all of the clusters come 
from the same underlying IMF (i.e., no scatter) the steep slope 
remains. The lack of correlations between mean IMF slope 
and cluster physical properties suggests that other potential 
sources of bias (e.g., mass segregation) contribute at a negli¬ 
gible level to our measurement of the mean. 

In contrast, it is more challenging to interpret the signif¬ 
icance of the upper bound on the scatter. The results of our 
inference suggest that a universal IMF is most probable. How¬ 
ever, there is a tail to the distribution that extends to > 0.1 dex, 
which exceeds the tail of the same distribution for the artificial 
clusters. We believe that this scatter is likely due to sources of 
uncertainty not included in our modeling that may artificially 
enhance the scatter relative to expectations from the artificial 
clusters. For example, mismatches between the stellar mod¬ 
els and real stars are not captured by the artificial cluster tests. 
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Table 3 

Select Literature High-Mass IMF Slopes 


Cluster 

(1) 

r 

(2) 

(M©) 

(3) 

LMC 

LH6 

L7±0.4 

85 

LH9 

L7±0.4 

55 

LH 10 

L7±0.4 

90 

LH38 

L7±0.4 

85 

LH 47/48 

L3±0.2 

50 

LH58 

L4±0.2 

50 

LH73 

L3±0.4 

65 

LH83 

L3±0.5 

50 

LH 114 

l.OiO.l 

50 

LH 117/118 

L6±0.2 

>120 

R136 

L3±0.1 

>120 

MW 

NGC 6823 

L3±0.4 

40 

NGC 6871 

0.9±0.4 

40 

NGC 6913 

Ll±0.6 

40 

Berkeley 86 

L7±0.4 

40 

NGC 7280 

L7±0.3 

65 

Cep OB5 

2.1±0.6 

30 

IC 1805 

L3±0.2 

100 

NGC 1893 

L6±0.3 

65 

NGC 2244 

0.8±0.3 

70 

NGC 6611 

0.7±0.2 

75 

Cyg OB2 

0.9±0.2 

110 

Tr 14/16 

L0±0.2 

>120 

h & X Per 

L3±0.2 

120 


Note. — High-mass IMF slopes from the massive star studies in the MW and LMC 
clusters as listed in Massey (1998) and Massey (2003). The dynamic mass range of the 
clusters extend from >1 Mq to the mass listed in column (3). 


However, the degree to which model-data mismatches con¬ 
tribute to the scatter is hard to quantify. 

Differential effects may also contribute to the observed 
scatter. For example, significant cluster-to-cluster variations 
in mass segregation or binarity could contribute to the scatter 
in ways not captured by the artificial clusters or our model. 
Similarly, we impose a single binary fraction for all masses, 
which may not be accurate, given that the most massive stars 
may be more likely to have equal mass companions than in¬ 
termediate mass stars (e.g., Sana et al. 2013; Kobulnicky et al. 
2014). At the same time, given the modest size of the upper 
bound, we expect these effects are, at most, fairly small. 

Overall, it is remarkable that the IMF slope is consistent 
with universal in this mass, metallicity, and cluster density 
regime. As summarized in Kennicutt (1998a), quantifying 
the degree of scatter is one of the most critical aspects of 
understanding the IMF, and ultimately how stars form out of 
dense gas. The intrinsic scatter we find suggests that individ¬ 
ual clusters should only minimally deviate from the fiducial 
value, and that extreme outliers (e.g., > 3-cr) are quite rare. 
The identification of such objects may represent insight into 
a mode of star-formation that truly differs from what we ob¬ 
serve in this regime. 

5.1. Comparison with Other High-Mass IMF Slopes in the 
Local Group 

It is important to place our results into context of other 
high-mass IMF studies in the Local Group (LG). Among the 
strongest evidence for the universality of the high-mass IMF 
in the LG comes from a series of papers that combine spec¬ 
troscopy and photometry of resolved, massive stars in clus¬ 
ters and OB associations in the LMC and MW (Massey et al. 
1995a,b; Massey 1998, 2003). 

To summarize, these studies found (a) no drastic differ- 



Figure 7. The number of stars predicted by various common IMF models, 
relative to Kroupa (dashed line). The vertical dotted lines indicate the mass 
range over which the high-mass IMF was measured in MSI. Our median IMF 
(solid red line) is essentially identical to that of Kennicutt (1983) (FKennicutt = 
1.5). However, we also find a small spread, which is indicated by the red 
band. Our IMF model indicates a modest deficit of stars >2 M© relative to 
the canonical Kroupa IMF. 


ences between the IMF slopes of very young, < 5 Myr star 
forming regions in the MW and LMC, despite large differ¬ 
ences in metallicity; and (b) that the resulting IMF slopes are 
broadly consistent with Salpeter/Kroupa. Combined, these 
studies represented a major milestone in our understanding 
of the high-mass IMF and its (in-)sensitivity to environment. 

Although re-analyzing these studies is beyond the scope 
of this paper, we can quantify what the ensemble high-mass 
IMF slopes are in the LMC and MW using data from Massey 
(2003). The purpose is to give a general sense of how different 
our M31 result is from other LG environments. To perform 
this analysis, we extracted the necessary data from Massey 
(2003) and have listed it in Table 3. We simply use the model 
from Equation 1 with = 1 for each object, i.e., we know 
they are all bone fide clusters, to infer the mean and variance 
for the MW and LMC clusters. We do not model possible age, 
mass, size dependancies. 

From this exercise we find mean values of Tuw = 1-16;^q}q 
and Flmc = L29;^o ii- t)oth cases, ay ^ 0.3-0.4. The large 
scatter in the posterior distributions is likely a reflection of 
underestimated uncertainties and/or systematics, as opposed 
to true physical variation. 

Massey (2003) emphasizes that the reported uncertainties 
are only statistical in nature. However, even these are likely 
under-estimated, as discussed in Weisz et al. (2013). Unfor¬ 
tunately, not all of these literature studies list both the number 
of stars and mass ranges needed to estimate degree of under¬ 
estimation of the random uncertainties. There are also issues 
of systematic uncertainties in stellar models that can signif¬ 
icantly change the masses of the most massive stars used in 
these studies (e.g., Massey 2011). 

However, even from this cursory exercise, we see that the 
MW, LMC, and M31 do not appear to have drastically dif¬ 
ferent IMF slopes. While this is quite remarkable given the 
diversity of environments, it should be tempered by the het¬ 
erogenous nature of this comparison. If high-mass IMF vari¬ 
ations are subtle, then homogenous and principled analyses 
are absolutely necessary to uncover them. At present, only 
our analysis in M31 is the result of a systematic study of the 
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Table 4 

Number and Fraction of Massive Stars for Different IMFs 


IMF 

m > 2Mq 

m > 8Mq 

V* > 20Mq 

V* > 50Mq 

> 90Mq 

Frac. >2 Mq 

Frac. >8 Mq 

Frac. >20 Mq 

Frac. >50 Mq 

Frac. >90 Mq 


10^ stars 

10^ stars 

10^ stars 

10^ stars 

10^ stars 







106 Mq 

106 Mq 

106 Mq 

106 Mq 

106 Mq 






(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

Kroupa 

6.9 

1.1 

3.0 

6.2 

6.3 

1 

1 

1 

1 

1 

Salpeter 

6.8 

1.0 

2.7 

5.4 

5.3 

0.97 

0.90 

0.87 

0.84 

0.83 

M31 

f- '7+0.1 
”•2-0.0 

”-°°-0.04 

'j+O.S 

^•^-0.2 

0 0 

bo 

Q Q+fi-9 
3.8 _o4 

0 90+0 04 

0 74+0.09 
”-/^-0.04 

n AA+011 
”•””-0.06 

r\ /crv+O.lS 
^•”±0.07 

n ^c/;+0.l4 
”•^”-0.06 

Kennicutt 

6.7 

0.82 

1.9 

3.5 

3.2 

0.87 

0.67 

0.57 

0.50 

0.46 

Miller & Scalo 

6.8 

0.66 

0.82 

0.81 

0.56 

0.85 

0.51 

0.23 

0.11 

0.08 

Scalo 

5.6 

0.20 

0.43 

0.72 

0.65 

0.54 

0.14 

0.10 

0.08 

0.07 


Note. — The number of expected stars above the specified mass limit per 10^ Mq. For these calculations, we have adopted a Kroupa (2001) IMF below 2 Mq, the specified IMF 
above 2Mq, and assumed lower and upper mass bounds for the total IMF to be 0.08 Mq and 100 Mq. Columns (7) - (11) give the fraction of stars formed above the listed mass, 
relative to the standard Kroupa IMF. 


high-mass IMF over a large number of young clusters. Any 
comparison of the high-mass IMF slopes of LG galaxies must 
come from similarly homogenous data and principled analysis 
in order to minimize the systematics that have been persistent 
in IMF studies for the better part of six decades. 

5.2. Practical Usage of High-Mass IMF 

Given the analysis presented in this paper and comparison 
with LG literature, we do not find strong evidence that a uni¬ 
versal Salpeter/Kroupa IMF is the best representation of the 
high-mass IMF slope in the LG. Although we only rule it 
out at a ^ 2 - cr level, there appears to be little basis for this 
canonical value being correct in the first place^^. Instead, for 
practical purposes the high-mass IMF can be represented by a 
single-sloped power-law with L45;^Q;Qg. Although this is only 
0.1-0.15 dex steeper than Salpeter/Kroupa it does make a dif¬ 
ference in the number of high-mass stars formed, and related 
quantities such as commonly used star formation rate (SFR) 
indicators. We discuss these implications further in §5.3. 

For general usage, we recommend adoption of the follow¬ 
ing forms of the the IMF over all stellar masses. With the M31 
high-mass IMF, a modified Chabrier IMF Chabrier (2003) has 
the form 


.. 0.086 -(log(m)-log(0.22))^ 

^(m) A(m) =- exp 2 x 0.572 

m 


am)=cm+™),r=1.45t!]l 

,for 1.0 < m < 100 Mq ,(3) 

and a modified Kroupa (2001) IMF has the form 

II 

p 

for 0.08 <m < 0.5 Mq 

^(m) = cm < F = 1.3 

for 0.5 < m < LOMq 

1 -p _ 1 4 C+ 0.03 

for 1.0<m< IOOMq. 


© • 

(4) 


Here, we have made two extrapolations to match conven¬ 
tional IMF definitions in the literature. First, we have extrap¬ 
olated the M31 IMF down to 1 Mq . For most practical pur¬ 
poses, this extension from 2 Mq to 1 Mq affects the resulting 
stellar mass spectra at the level of a few percent. 

Second, we have extrapolated the M31 IMF up to 100 Mq. 
While our data do not provide good constraints above ^ 25 

Salpeter (1955) has had remarkable staying power, despite the fact that 
the power-law index of F = 1.35 was derived assuming the MW was 6 Gyr 
old. A more realistic 13 Gyr age of the MW yields a value of F = 1.05, which 
is considered an extreme outlier by modern standards (e.g., Salpeter 2005; 
Zinnecker 2005). 


Mq, there is little other information to go on for mass spec¬ 
trum of the most massive stars. Thus, as is commonplace in 
the literature, we suggest extrapolating the M31 IMF up to the 
highest masses, but recognize that systematic spectroscopic 
searches for the massive stars are absolutely critical for un¬ 
derstanding the IMF of the most massive stars (e.g., Massey 
et al. 1995b; Massey 2003). 

Finally, we have assumed that the sub-solar IMF in M31 
is similar to that of the Galaxy. Although not confirmed by 
direct star counts, Conroy & van Dokkum (2012) have shown 
that spectral features sensitive to the low-mass IMF yield a 
result that is similar to the low-mass Galactic IMF. 

5.3. Broader Implications 

In Figure 7 we illustrate the mass spectra of high-mass stars 
predicted by various commonly used IMF models. The plot 
shows the number of stars expected relative to a Kroupa IMF 
(F = 1.3) as a function of mass. For fair comparison, the 
slopes have all been normalized at 2 Mq, the lower main se¬ 
quence mass limit of the PHAT data, and extended to 100 Mq 
on the upper end. We chose this normalization, as opposed to 

1 Mq to discuss out ability to distinguish different literature 
forms of the IMF over the dynamic mass range of the data, 
which is indicates by the vertical dashed lines. 

The M31 high-mass IMF predicts the formation of fewer 
stars than Kroupa at all masses. The red solid line presents 
our median IMF model and shows that the fractional deficit 
varies from ^0.9 at ^ 10 Mq to ^ 0.7 at 100 Mq. The 
red shaded region reflects the 68 % confidence interval in the 
number of predicted high-mass stars due to uncertainty on the 
mean high-mass IMF slope in M3L At most masses, the un¬ 
certainty deviates from the median by ±^5%, but increases 
to ±15% for stars > 50 Mq. 

The M31 high-mass IMF is similar to other high-mass IMF 
models in the literature. It is particularly close to the Ken- 
nicutt IMF (F = 1.5; Kennicutt 1983). It is also the same as 
the IMF of Miller & Scalo (1979) over the 1-10 Mq range. 
The steepness of the Miller & Scalo (1979) IMF above 10 
Mq is known to be too extreme (e.g., Kennicutt 1983; Ken¬ 
nicutt et al. 1994). Further, the high-mass IMF from Scalo 
(1986) has a similar slope (F = 1.6) above ^ 4 Mq, but a 
much steeper slope below it. Finally, the M31 high-mass IMF 
predicts fewer high-mass stars than a Salpter IMF for masses 

>2Mq. 

We further quantify differences in these IMF models in Ta¬ 
ble 4. Here, we have computed the expected number of stars 
per 10^ Mq formed at selected stellar masses assuming a 
Kroupa (2001) IMF below 2 Mq and the listed IMFs above 

2 Mq . This table solidifies many of the above points regard- 
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Table 5 

Updated Broadband SFR Indicators 


Band 

Lx Units 

log(CU 

log(G) 

M*(M31)/M*(Kroupa) 



(Kroupa) 

(M31) 


(1) 

(2) 

(3) 

(4) 

(5) 

EUV 

erg s“^ {vL^) 

43.29 

43.18!g:S5 

1 OQ“*“t).06 

A-^0-0.12 

NUV 

erg s“^ {vLj,) 

43.14 

43.04!S;S3 

1 24+0-0^ 

Ha 

erg s“^ 

40.91 

40.73!S:S^ 

1 52“''®'^^ 


Note. — A comparison of commonly used SFR indicators for a Kroupa and M31 IMF 
values. The values of L^, Cx, and M* follow Kennicutt & Evans (2012) and Equation 5. 


ing differences in commonly used high-mass IMF models. It 
also highlights an important implication for core collapse su- 
pemovae, whose progenitors are believed to have masses > 
8 M 0 (e.g., Smartt 2009). For example, the median M31 high- 
mass IMF model predicts ^25% fewer stars with masses > 
8 M 0 compared to a Kroupa IMF. 

Finally, we consider the implications for commonly used 
SFR indicators (e.g., Kennicutt 1998b; Kennicutt & Evans 
2012). Using the Flexible Stellar Population Synthesis code 
(FSPS; Conroy et al. 2009, 2010), we compute the luminosity 
to SFR conversion coefficients for GALEX far- and near-UV 
luminosities, along with Hcf, which we list in Table 5. To do 
this, we have assumed a constant SFH over the last 1 Gyr, the 
indicated IMF, mass limits from 0.08 to 100 M 0 , and case 
B recombination to convert the number of ionizing photons to 
Ha luminosity. From Kennicutt & Evans (2012) the canonical 
conversion from luminosity to SFR can be written as 

logM*(M©yr“b = logL;^ - logCx (5) 

where is the SFR, is the luminosity over the bandpass 
X, and Cx is the conversion constant over the same wavelength 
range. We have provided updated values for Cx in Table 5. 

From column (5) of Table 5, we see modest differences in 
the SFR indicators due to an M31 high-mass IMF. The varia¬ 
tions are such that for a fixed luminosity, using the M31 high- 
mass IMF results in SFRs that are a factor ^1.3-1.5 larger. 
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